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ABSTRACT 

Recent Chandra and Newton observations indicate that there are two-temperature 
components (T ~ 8 keV, 0.8 keV) of the diffuse x-rays emitted from deep inside 
the center of Milky Way. We show that this can be explained by the existence of 
sterile neutrinos, which decay to emit photons that can be bound-free absorbed by 
the isothermal hot gas particles in the center of Milky Way. This model can account for 
the two-temperature components naturally as well as the energy needed to maintain 
the ~ 8 keV temperature in the hot gas. The predicted sterile neutrino mass is between 
16-18 keV. 
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1 INTRODUCTION 

Recently, a large amount of diffuse x-ray data have been 
obtained by Chandra, BeppoSAX, Suzaku and XMM- 
Newton, giving a complex picture near the Milky Way center 
jMuno et al.ll2004[Ro~ckefeller et alj20q4ISenda et al.ll2002l; 
Hamaguchi et all 120071 ; ISidoli et alj Il999l ; ISakano etaU 
20041 ). The data indicates that there exists a high tem- 



perature (~ keV) hot gas near the Milky Way center. 



iKaneda et ail l|l997l) suggested that the existence of two- 
temperature component s of the hot ga s can e xplain the ob- 
served x-ray spectrum. ISakano et al.l (|2004T ) analysed the 
XMM-Newton data to get 1 keV and 4 keV hot gas com- 
ponents in Sgr A East, a supe rnova remnant loca ted close 
to the Milky Way center. Later iMuno et al.1 (|2004l ) used the 
data from Chandra to model the temperature of the two 
components within 20 pc as 0.8 keV and 8 keV. The tem- 
perature of the soft component (Ti ~ 0.8 keV) can be ex- 
plained by 1 percent of kinetic energy by one supernova 
explosion in every 3000 years near the Milky Way center, 
corre sponding to power ~ 3 x 10 36 erg s _1 |Muno et al.l 
120041 ). However, the high temperature of the hard compo- 
nent (Ta ~ 8 keV) cannot be explained satisfactorily. The 
power needed to sustain the high temperature is about 10 40 
erg s _1 . Moreover, the emission of the hard component is dis- 
tributed much more uniformly than the soft component and 
the intensity of the two components are correlated, which 
suggest that they a re produced by related physical processes 
IjMuno et al.ll2004T ). In this article, we present a model to ex- 
plain the two-temperature composition, as well as the high 
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temperature and uniform emission of the hard component 
obtained by Chandra within 20 pc. We assume that there 
are sterile neutrinos with rest mass m s ^16 keV in the deep 
galactic center and their decay photons continuously supply 
the energy of the hard x-ray emission. The sterile neutrino 
halo is located at the center of a 20-pc radius gas cloud 
which is heated by the decayed photons. Most photons are 
absorbed near the center (~ — 1 pc) of Milky Way and 
the energy is subsequently transferred to the surrounding 
gas. The major heating mechanism is the bound-free col- 
lisions between the decay photons and the ions in the hot 
gas. In this optically thick region, the hot gas particles are in 
photoionization equilibrium. The high metallicity of the gas, 
Zmetai ^2 — 3 solar metallicity, enhances the heating rate 
of the gas and provides enough energy to sustain the high 
temperature of the hard component. The energy absorbed 
in the central region is transferred to the outside optically 
thin region (1 — 20 pc) by collisions among the electrons 
to share their energy (collisional equilibrium). Also, we as- 
sume that the soft and hard components are in equilibrium 
with different uniform temperatures and they are bounded 
hydrostatically. 

Although the recent MiniBooNE data challenges the 
LSND result th at suggests the existence of eV scale 
sterile neutrinos (|Aguilar-Arevalo et al.l 120071 ) , more mas- 
sive sterile neutrinos (eg. keV) may still exist. The fact 
that active neutrinos have rest mass implies that right- 
handed neutrinos should exist which may indeed be mas- 
sive sterile neutrinos. The existence of the sterile neutri- 
nos has been i nvoked to explain many p henomena such 
as reionization llHansen and Haimanl [2004). missing mass 
|Dodelson and Widrow 1994h and the high temperature of 
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the hot gas in clusters l|Chan and Chu 2007). Therefore, it 
is worthwhile to discuss the consequences of the existence of 
massive sterile neutrinos, which may decay into light neu- 
trinos and photons. 

The existence of the small size keV ster il e neu - 
trino halo is first suggested by IViollier et al.l (| 1993h . 
The size of a self-gravitating degenerate sterile neu- 
trino halo depends on m 3 and total mass M 3 : R a = 
O.OOO6(Af s /lO 6 M )" 1/3 (m s /16 keV)" 8/3 pc. Including the 
contribution of the baryons, the size will be even smaller. 
The size of the ste rile neutrino halo is upper bounded by 
R s < 0.0005 pc (|Schodel et al.l |2002| ). This size is very 
small compared to that of a galaxy and therefore the ster- 
ile neutrino halo will hide dee ply inside the galactic center 
l|Munvaneza and Vio llicr 2002). Sterile neutrinos may decay 
into active neutrinos and become a strong energy source to 
galaxies and clusters. The decays of keV order sterile neu- 
trinos m ay also help to solve the cooling flow problem in 
clusters l|Chan and Chull2007l ). 



2 BOUND-FREE ABSORPTION MODEL 

The sterile neutrinos at the center will decay into active 
neutrinos with decay rate F by the following process: 



v a + 7- 



(1) 



We assume that the energy of decay photons E a « m a /2 
is greater than 8 keV because the energy of each photon 
must be greater than the energy of each electron in the hot 
gas in order for the l atter to gain energy from the photons 
l|Chan and Chu|[2007l ). The decayed photons come from the 
volume omission of the entire sterile neutrino halo. The dis- 
tribution of photon energy has a characteristic width deter- 
mined by the Fermi momentum of the sterile neutrinos pf, 
which is very small compared with the rest mass of the sterile 



neutrinos, pf/t 



10" 



Therefore, we can approximate 



the photon spectrum as monochromatic with energy E a . 

Since we have not detected any strong lines of such 
high energy photons from the Milky Way center, the optical 
depth for decayed photons must be much greater than 1. 
Therefore, the total energy emitted by the decayed photons 
must be equal to the total energy gained by the electrons 
(~ 10 40 erg s- 1 ): 



NsT(E a - E,)P, ^ 10 40 erg s" 



(2) 



where N a is the total number of sterile neutrinos, Ei and 
Pi are the ionization potential and probability of photon 
absorption by i th type ions in the hot gas: 

0,i<Jbf,i 



P 



J2i a i a bf,, 



(3) 



where is the ratio of the number of i* type ions to the 
total number of ions at 0.8 keV temperature (the number 
density of the soft component is about ten times more than 
that of the hard component if they are in equilibrium). The 
absorption cross section of the i th type i ons is largest for H- 
like a nd He-like ions, which is given by |Daltabuit and Coxl 
[19Z2): 



(Tbf,i 



10 ath,- 



Eth.i 



where (Tth,ii Q ii Eth,i and s, are fitted parameters of a 
particular i th type ion. The total effective cross section is 
y\ a,i(Jbf,i- The metallicity Z mota i near the Milky Way cen- 
ter ca n reach 2 — 3Z me taX,@ (solar metallicity) (ISakano et al.l 
2004). The metal ions with the largest cross section include 
Si (Zsi = 8.9Z Si , Q ), S (Z s = 2.7Z s>e ), At (Z Ai = 1.8Z Ar , ), 
Ca (Z C a = 2.5Z Ca .,e ) and Fe (Z Fc = 3.8Z Fo , ) (see table 1) 
l|Sakano et al.1120041 ). Assuming Z meta i = 3-Z me tai,0 for other 
metals and V = (3 — 6) X 10~ 19 s -1 ( the best fit value tha t 
can solve the cooling flow problem) (|Chan and Chull2007l ). 
we can obtain N 3 and the total mass of the sterile neutrino 
halo (M s = N s m B ) by Eq. (2) for different E s (see table 
2). The value of M a obtained (~ 1O 5 M ) is far below the 
upper bound of the total mass nea r the Milky Way center 
(2.6 x 10 6 Af Q ) jSchodel et alj|2002ft . 

To constrain the value of E a , we can consider the upper 
limit of x-ray luminosity in the 2 — 10 keV b and near the 
Milky Way center, which is 3 x 10 36 erg s _1 l|Muno et al.1 
l2004h . We can obtain the lower limit of the optical depth 
Tiowcr by: 

N a FE a e' T «S 3 x 10 36 erg s _1 , (5) 

where r is the optical depth for decayed photons which is 
given by: 



t= I n(r) ^ aj(Tbf,idr. 



(6) 



A higher energy of the decayed photons will result in a 
smaller cross section and optical depth. The mass density 
profile of Milky Way center can be modeled by a power law 
p ~ r" 1 ' 8 l|Schodel et alj|2002h . For simplicity, we assume 
that the number density of hot gas near the center has an 
isothermal profile: 

nW=no( I ^) 2 . (7) 

The average number density within 20 pc of the center is 
about 1 cm -3 which corresponds to no = 133 cm -3 . In table 
2, we can see that r < ri owcr for E s ^ 10 keV. Therefore, 
only E 3 = 8 — 9 keV can be possible for decaying sterile 
neutrinos as the heating source (see Fig. (1)). 

Given this power input (10 40 erg s^ 1 ), we can esti- 
mate the equilibrium temperature of the gas at the center 
of Milky Way. The ionizing photons are continuously sup- 
plied and the heating rate is quite constant as the decay of 
sterile neutrinos is a slow process (half-life of order Hubble 
time). The energy absorbed will then be transferred to the 
other gas particles mainly by conduction. The power loss by 
Bremsstrahlung radiation of a hot gas with temperature T 
is A = 1.4 x 10 27 riiT 1/2 erg s -1 cm -3 , and the power by its 
adiabatic expansion is 



W = PV 2/3 c s = PV 



r 2/3 



(lkT\ 
V m a J 



1/2 



(8) 



■(41 



where P, V, c s , 7 and m g are pressure, volume, sound speed, 
adiabatic index of the hot gas and mean mass of the gas par- 
ticles respectively. The energy loss due to Bremsstrahlung 
is negligible, being only 0.3 percent of the adiabatic cooling 
(M uno et al. 2004). We can therefore solve for the equilib- 
rium temperature of the gas using Eq. (8). The temperature 



maintained in this process is given by 
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Table 1. Metal abundances we have used in the model. 



Element 


Atomic number Z 


Metallicity Z me tal 


ai(10- 4 )(H-like) 


ai(10- 4 )(He-like) 


C 


6 


3 


6.9 


2.7 


N 


7 


3 


1.7 


0.84 





8 


3 


16 


9.6 


Ne 


10 


3 


2.5 


1.9 


Mg 


12 


3 


1.0 


0.89 


Si 


11 


8.9 


.3.0 


2.8 


S 


16 


2.7 


0.43 


0.42 


Ar 


18 


1.8 


0.066 


0.066 


Ca 


20 


2.5 


0.056 


0.056 


Fe 


26 


3.8 


1.8 


1.8 




Figure 1. The allowed N B and m s by using the constraints in Eq. (2) (dashed line), Eq. (5) (solid line) and by the temperature of the 
hard component (E B ^ 8 keV) (dotted line). 



T « V 



-4/9 



w 

n t k 



2/3 



yk 



1/3 



(9) 



If W = 10 40 erg s" 1 within 20 pc, then T « 8 keV. In 
fact, the situation is similar to energy absorption in x-ray 
emitting nebulae. However, in a normal x-ray nebula, the 
x-ray source is usually far away from the gas cloud, whereas 
in our model the source is located at the center of the gas 
cloud. Furthermore, the energy absorption in a normal x-ray 
nebula is about 10 30 erg s -1 so that the tem perature of the 
nebula is about eV order (|Leahv et al.l ll994T). In our model, 
the energy absorption is 10 40 erg s _1 , which can maintain 
the temperature of the hard component at 8 keV. 



3 TWO-TEMPERATURE COMPONENTS 

When the gas particles are in thermal equilibrium, 

f(n g ,T) = A(n g ,T), (10) 

where n g is the number density of the gas particles, T and 
A are the heating and cooling rates respectively. It is pos- 
sible to have multi-temperature components in the hot gas, 
because there may be more than one solutions (n g ,T) sat- 
isfying Eq. (10). The criteri on for stability of the solution is 
ijBowers and Deemind[l984l ): 



(11) 



where p and P are the mass density and pressure of the hot 
gas. Since the relaxation time of electron-electron collisions 
is less than that of ion-electron collisions, the temperature of 
electrons T e may be different from the temperature of ions 
Ti in general. The total heating rate inside the gas cloud 
(r ^ 7? m 20 pc) is given by: 



(12) 



where To is a constant that depends on the size of the 
region. The cooling rate of x-ray emitting hot gas by 
Bremsstrahlung radiation is: 



A ~ A(,n e n g CbiZ^e 



-2I i /kT i rpl/2 



(13) 



where Ao is a constant, Ii and Zi are the ionization en- 
ergy and charge of the i th type ions respectively. Since 
the entire gas cloud is optically thin for r ~ pc and 
heavy metal i ons a re concentrated at the central region 
l|Sakano et alj [2004), the bound- free absorption of the 
broad-band Bremsstrahlung photons is suppressed (r ^ 1) 
in a large fraction of the volume of the cloud, so that most of 
the Bremsstrahlung photons can escape the gas cloud eas- 
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Table 2. The values of N s , M s and r for E a 



= 8-14 keV. 



E s (keV) 


7V S (10 66 ) 


M S (1O 5 M ) 


J2i ^o-^lO" 23 cm 2 ) 


T 


^lowcr 


8 


13 


1.9 


1.77 


14.6 


10.0 


9 


11 


1.7 


1.26 


10.4 


9.93 


10 


9.1 


1.(5 


0.921 


7.61 


9.87 


11 


1.6 


0.9 


0.685 


5.74 


9.30 


12 


3.6 


0.76 


0.536 


4.43 


9.12 


13 


2.9 


0.67 


0.422 


3.48 


8.90 


14 


2.4 


0.61 


0.338 


2.79 


8.90 



ily. When the hot gas particles are in thermal equilibrium 
r = A, and the pressure is given by: 

rpl/2 v-^ 

J k e vLi a ' ff ^' 



P = Po-, 



(14) 



where P = (1 keV) 1/2 fcf /A and T kcV = T e /(1 keV). 

Fig. (2) shows the relation bet ween P and T e by using 
the mean metallicity within 20 pc (|Muno et al.ll2004h . We 
notice that for some values of P there are two different values 
of T e . By putting the heating and cooling rates into Eq. (11), 
we get the criterion for stable solutions: 



J2i a i°bf,i 



5>Z|e^ (| -ft) 



dai 



where f3i = 2U /kTi. We impose a free parameter \x such 
that T e = fj,Ti. We define A and B to be the expressions on 
the left and right hand sides of the inequality in Eq. (15) 
respectively. In Fig. (3), we plot A — B against T e for two 
values of fi. When fj, = 3.2 and fi = 6.9, 0.8 keV < T e < 
4 keV and 1.8 keV < T e < 8 keV are unstable solutions 
respectively as A — B > 0. We also indicate the unstable 
regions in Fig. (3), and we notice that a two-temperature 
phase may exist with T e ^ 8 keV and T e ^ 0.8 keV for 
3.2 ^ /i ^ 6.9. In Fig. (2), we can see that if there are less 
large- Z metal ions (especially iron) in the hot gas, more hot 
gas particles will shift to the lower temperature phase and 
only one equilibrium solution may be obtained. Clearly, our 
model can account for the two-temperature structure of the 
hot gas near the Milky Way center. 



4 DISCUSSION 



To ex plain the origin of the hard component, iMuno et al.l 
i|2004h make use of magnetic reconnection driven by the tur- 
bulence that supernovae generate in the interstellar medium. 
Magnetic reconnection can heat the hot gas to kT ~ 

Bcentei /8lVn s . Foi Tig ~ 0.1 Cm" 3 , -Bcontcr ~ 0.2 HlG, 

kT ~ 8 keV. However there is not enough evidence to sup- 
port whether this mechanism ca n maintain the h igh tem- 
perature of the hard component (|Muno et al.ll2004l ). 

In our model, we have assumed that there exists a ster- 
ile neutrino halo with m s = 16 — 18 keV in the Milky Way 
center, which decay to emit 7 with life-time of cosmological 
order. It provides a large amount of energy to the hot gas 
and maintains the extremely high temperature. The bound- 
free collisions provide enough energy to the two different 
temperature components and maintain their temperatures. 
At the same time, a stable two-temperature structure in the 



hot gas can be explained by this heating mechanism natu- 
rally. The uniform emission of the soft and hard components 
suggests that they may come from similar physical processes 
l|Muno et al.ll2004l ). In our model, both components indeed 
share the same source of energy - the 8-9 keV photons emit- 
ted by the decays of sterile neutrinos. In our model, the 
sterile neutrinos may not be the major component of dark 
matter. Therefore, any bounds on m s assuming they are the 
major dark matter candidate does not constrain our model 
severely. The heating rate in the Milky Way center is time 
dependent as there is a decreasing number of sterile neutri- 
nos. Therefore, if two galaxies have similar chemical compo- 
sitions, the heating rate is greater for the large redshift one, 
has more particles in the higher temperature compo- 
nent. We therefore predict that the hard component of the 
x-rays would be stronger for large redshift and metal-rich 
galaxies. Moreover, if a galaxy has lower metallicity, then 
only a single temperature component may be observed in- 
stead of two (see Fig. (2)). 
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Figure 2. Pressure (P) versus temperature (T) of the hot gas, given by Eq. (14) for Zp e = 0.71 (solid lines) and Zp e = 0.14 (dashed 
lines) for two different values of fi. 
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Figure 3. The difference between the left hand side (A) and right hand side (B) of Eq. (15) versus T e for fj, = 2.7 (dashed line), \x = 3.2 
(solid line) and fi = 6.9 (dotted line). 
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